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ABSTRACT 

We present a numerical investigation of the coronal evolution of a coronal mass ejection (CME) 
on 2005 August 22 using a 3-D thermodynamics magnet ohydro dynamic model, the SWMF. The 
source region of the eruption was anemone active region (AR) 10798, which emerged inside a coronal 
hole. We validate our modeled corona by producing synthetic extreme ultraviolet (EUV) images, 
which we compare to FIT images. We initiate the CME with an out-of-equilibrium flux rope with an 
orientation and chirality chosen in agreement with observations of a Ha filament. During the eruption, 
one footpoint of the flux rope reconnects with streamer magnetic field lines and with open field lines 
from the adjacent coronal hole. It yields an eruption which has a mix of closed and open twisted 
field lines due to interchange reconnection and only one footpoint line-tied to the source region. Even 
with the large-scale reconnection, we find no evidence of strong rotation of the CME as it propagates. 
We study the CME deflection and find that the effect of the Lorentz force is a deflection of the 
CME by about 3° / R© towards the East during the first 30 minutes of the propagation. We also 
produce coronagraphic and EUV images of the CME, which we compare with real images, identifying 
a dimming region associated with the reconnection process. We discuss the implication of our results 
for the arrival at Earth of CMEs originating from the limb and for models to explain the presence of 
open field lines in magnetic clouds. 

Subject headings: Sun: corona — Sun: coronal mass ejections (CMEs) — magnetohydrodynamics: 
MHD — Sun: magnetic topology 



1. INTRODUCTION 

\ Coronal mass ejections (CMEs) are one of the lead- 
ing causes of space weather, especially when they have 
organized southward directed magnetic field. Therefore, 
it is of great importance for space weather forecasting 
to understand how their direction of propagation and 
their orientation at 1 AU relate to properties on the 
solar disk (location of source region, orientation of the 
polarity inversion line, etc.). The deflection of CMEs 
in the latitudinal direction has been observed and re- 

] ported since the launch of the first sp ace corona^raphs 
in the 1970s and 1980s. For example, MacOueen et al.l 
(p. 986) reported an average deflection of 2.2° towards 
the equator for 29 CMEs during solar minimum (1973- 
1974), while they found no systematic deflection for 19 
CMEs during solar maximum (1980). With the launch of 
the Solar- Terrestrial Relations Observatory (STEREO) 
in 2006, there has been a renewed focus on the deflec- 
tion of CMEs in the corona and the heliosphere, thanks 
to stereo scopic measurements from the two spacecraft 
(e.g. see: iKilpua et al.ll2008l : iLiu et al.l[20Tol : iBvrne et al.l 



The presence of coronal holes is known to affect the 
direction of propagation of CMEs (|Plunkett et al.|[2QQTI : 
[Gopalswamy et al. 2009). Howe ver, the exact cau se of 
this influence is uncertain. ICremades et al.l ()2006f ) and 
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iGopalswamv et ahl ()2009f ) have used an ad hoc, fictitious 
force exercised by coronal holes on CMEs, which is di- 
rectly dependent on the coronal hole area and the dis- 
tance of the CME source region from the coronal hole. 
This force might be due to a ^^pushi ng" of the CME b y 
the fast wind f r om th e coronal hole (jWang et al.ll2004f ). 
iPlunkett et al.l (|2001f ) proposed that the effect of coro- 
nal ho les is due to strong magnetic fields, iFilippov et al.l 
(|2001f ) proposed that the non-radial motion is due to 
the g:uidings: action of the co ronal m agnetic field, while 
lAulanier et al.l (j2010D and IShen et al.l (|201lf ) recently 
proposed that magnetic pressure gradient acting on a 
CME results in a net force directed along the gradient of 
the magnetic pressure. Obviously, magnetic forces, such 
as the Lorentz force could also cause the coronal hole 
to exercise a force on the CME. Magnetic forces are ex- 
pected to affect strongly the evolution and propagation 
of a CME in the corona while the effect of the fast so- 
lar wind from the coronal hole would continue into the 
heliosphere and be more gradual. 

Similar effects may result in the rotation of an erup- 
tion. One of the leading causes of CME rotation is the 
kink instability (jTorok fc KliemI l2003f ) but recent nu- 
merical efforts have focused on other sources of rota- 
tion, such as reconnectio n with the background magnetic 
field (Cohen et al."2010'), the effect of the Lorentz force 
(Isenberg & Forbes 2007; Shiota et al. 2010) and its asso- 
ciation with shearing motions (Lynch et al. 2009). These 
numerical works point towards a near universal rotation 
for CMEs not initially aligned with the heliospheric cur- 
rent sheet (HCS). While statistically speaking, CMEs 
rotate towards the HCS (Yurchyshyn 2008), it appears 
from observations that some CMEs do not rotate, even 
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though they are ori^inany not aliened with the HC S (see 
some cases m lYurchvshvnll2QQ8l : lYurchvshvn et alJ [2QQ9) . 

Eruptions from active regions which are inside a coro- 
nal hole present a perfect configuration to study the ef- 
fect of coronal holes on CMEs. F or example, the ficti- 
tious force of Cremades et al.l ()2QQ6 ) is inversely propor- 
tional to the square of the distance between the source 
region and the coronal hole, and it results in a infi- 
nite force for a CME originating from such an active 
region! Active regions can form inside coronal holes 
and, then, the y oft en develop into anemone active re- 
gions (gtubMa^elalJ [lE^ Anemone 
active regions were origi nally observed in sof t X-ray and 
named X-ray fountains (jSheelev et al.l [19751 ) due to the 
structure looking like a fountain with loops emerging in 
all directions from the center. Almost identical struc- 
tures are observed in extreme ultraviolet (EUV) and 
in chromospheric lines (Shibata et al. 2Q0Z). In a sur- 
vey of active reg ions observed by th e soft X-ray tele- 
scope (SXT, see: iTsuneta et aLlflDQlD o nboard Yohkoh 
(jOgawara et al.lll991f ). lAsai et al.l (j2QQ8f ) found that as 
many as 25% of active regions observed in 1991-1992 
appeared as anemone active regions at one time during 
their evolution, and that there is a near equivalence be- 
tween an active region having an anemone structure and 
being inside a coronal hole. 

Anemone active regions have been primarily studied as 
the source of X-ray jets since the emergence of a bipo- 
lar active region inside unipolar open magnetic fields 
naturally yield such phenomena (jShibata et al.l [19941; 
IVourlidas et al.l [19961 ). Of greater importance for space 
weather are the instances of CME s and filament erup- 
tions from anemone active regions (iChertok et al.lIM)! 
lAsai et alj 120091 : iBaker et al.l [20091). iLiul (|2QQr) found 
that eruptions from active regions inside open magnetic 
field (which contain anemone active regions) are statisti- 
cally faster on average than other eruptions (even origi- 
nating from under the heliospheric current sheet). Most 
authors have explained this result by invoking two facts. 
Since the eruption originates from a region of low-density, 
fast solar wind, the fast wind can "push" the CMEs to 
a faster speed. Closed field in the low corona can hinder 
the eruption of the CME be cause of th e Lorentz force 
(see for example the model of lChen|[l996[ ) and in the ab- 
sence of closed field lines, the CME can erupt without 
being strongly decelerated. 

In the present work, we investigate the effect of the par- 
ticular magnetic structure of anemone active regions on 
the evolution and coronal propagation of a CME. We fo- 
cus on the first of two eruptions in 2005 August 22 from 
active region NOAA 10798. This CME was associated 
with the eruption of a southward-directed Ha filament 
and reached a speed of about 1200 k ni s~-^ (for a com- 
plete overview of the observations, see lAsai et al.|[2009 ). 
The second eruption was faster (~ 2400 km s~-^) and it 
occurred about 16 hours after the first one from the same 
active region. The eruption of a preceding CME from 
the same active region is expected to signific antly change 
the m agnetic topology of the source region (jLugaz et al.l 
l2010f ) and also modify the backgrou nd solar wind into 
which the second CME propagates (jLugaz et al.l l2005l 
[2007), which would complicate a numerical investigation 
of the second CME. For this reason, we focus on the first 



CME. Our investigation is based on a numerical magne- 
tohydrodynamic (MHD) model par t of the spac e weather 
modeling framework (SWMF, see: iToth et all 12005). In 
order to validate our simulation with EUV observations 
and to reproduce more accurately the lower part of the 
corona, we use the low corona (LC) model recently devel- 
oped bv lDowns et al.l ()2010i) . which captures the energy 
balance of the coronal plasma. 

The organization of the paper is as follows: in Sec- 
tion [21 we briefiy summarize the main features of the LC 
and flux rope models important for this study. Then, 
we discuss the magnetic topology of the anemone active 
region before presenting a comparison of the pre-event 
corona as observed in EUV with the modeled corona in 
synthetic EUV. In Section [H we follow the initial phases 
of the eruption and explain its interaction with the adja- 
cent magnetic flux systems. In this section, we also dis- 
cuss the CME aspect in real and synthetic coronagraphic 
images. We study the CME rotation, deformation and 
deflection in Section [H We discuss our results and con- 
clude in Section m 

2. NUMERICAL MODELS AND PRE-EVENT MODELING 
2.1. Low Corona Model 

The simulation is done using the Low Corona (LC) 
component of the SWMF. The LC model includes ra- 
diative losses, fleld-aligned electron heat condition and 
empirical heating tern is in the energy equation as de- 
tailed in Do wns et aJI ()2010 ) . It has been recently used 
to study the nature of the extreme ultraviolet (EUV) 
wave obse rved by the two S TEREO spacecraft in 2008 
March 25 (|Downs et al.ll20lT[ ). The initial magnetic fleld 
and boundary conditions are provided by a flnite differ- 
ence solu tion of the potential fie ld source surface method 
fPFSSM: [Altschuler et al.lll977[ ) for the synoptic magne- 
togram of Carrington ro tation 2033 as ob served by the 
SOHO/Up i instrument (jToth et al.ll2011[ ). As noted in 
lAsai et al.l (|200 9). AR 10798 was rapidly evolving as it 
crossed the Sun central meridian and therefore, the syn- 
optic map does not fully capture the complexity of the 
active region as revealed by observations on August 21 
and 22. However, on these two days, the AR was al- 
ready too close to the eastern limb of the Sun, so that 
the daily line-of-sight magnetogram (shown in the left 
panel of Figure 1) cannot be adequately used to con- 
strain our model. As discussed below, even using the 
synoptic map, the MHD model successfully reproduces 
the important featu res of the anemone a ctive region. 

As described in iDowns eTall ([201(1. we followed 
the previous works of Lionello et al. (2001) and 
iLionello et al.l (2009) to widen the transition region by 
modifying the ratio of heat conduction to radiative cool- 
ing at temperatures below 300,000 K. The heating model 
and boun dary conditions are the same as the ones de- 
scribed in iDowns et al.l (|2011[ ). i.e. we use the chromo- 
spheric boundary conditions. The grid is spherical with 
non-uniform radial scales near the transition region (min- 
imum uniform dr ~ 230 km) with a resolution on the 
solar surface of 1.4° x 1.4°. In addition, we further refine 
around AR 10798 by 2 levels of refinement as follows. 
The solar surface is refined twice in a box of about 20° 
latitude and 40° longitude around the active region over 
a height of 0.4 Rq and the region around the top of the 
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anemone structure (null point) is refined by an additional 
factor of 2 in a twice smaller box centered at 1.15 R©, 
resulting in a smallest angular resolution of 0.35° x 0.35°. 
The grid around AR 10798 is shown in the top left panel 
of Figure 2. 

This version of the LC model includes a uniform 
Parker-like solar wind model but not a bimodal solar 
wind model as in th e mod els of iCohen et al.l ()2007[ ) and 
Ivan der Hoist et"al] ()2010l ). Here, the solar wind accel- 
erates and becomes faster than the fast magnetosonic 
speed between 6 and 14 R© and reaches a value of 300- 
350 km s~^ at the outer boundary of the domain (24 R©). 
While it is not realistic for the fast wind originating from 
coronal holes, it is the approximate speed and accelera- 
tion profile expected for the slow solar wind. We only 
focus on the coronal evolution of the CME, therefore re- 
solving a bimodal solar wind is not critical to our study. 
Since the solar wind is uniform, we can expect a neg- 
ligible amount of deflection and rotation due to hydro- 
dynamical effect (velocity shear, non- uniform drag). It 
allows us to focus solely on the effect of magnetic forces 
and reconnection on the CME evolution. Compared to 
previous studies of the evolution of CMEs in the corona 
with the SWMF (IIjiga^^laD^^ ICohen et al.ll2010l : 
lEvans et al .1120 111 ), using the LC model has a number of 
advantages. The inclusion of the thermodynamics makes 
the treatment of the lower corona more realistic. Us- 
ing the chromospheric boundary conditions, the density 
and temperature of the plasma are based solely on the 
magnetic topology and the part icular choice o f hea tings; 
function. Also, as discussed in iDowns et all (j2010f ). it 
results in a less potential steady-state solution for the 
background coronal magnetic field. Last, it gives us the 
opportunity to validate our model with a direct compar- 
ison to extreme ultraviolet (EUV) observations. 

2.2. CME Model 

To model the CME, we employ a semi-circular flux 
rope presc ribed by a given total toroidal current , as in the 
models bv iTitov fc Demoulin (1999) and Rous sev et al.l 
(|2003f ). An azimuthal current is also added at the sur- 
face of the flux rope in order to construct a force-free 
magnetic field inside the flux rope. The toroidal current 
is largely dominant, creating a very highly twisted flux 
rope. A more complete description of this implementa- 
tion of the flux rope model ca n be found in Lugaz et al. 
120 07) andlEvans et al.l (|2QTiI) . Because in this version of 
the iTitov fc DemoulinI ()1999f ) flux rope, there is neither 
sub-surface line current nor sub-surface magnetic charges 
to generate strapping fleld for the flux rope in the corona, 
the flux rope is not expected to rotate as found for exam- 
ple in Isenberg fc Forbes (2007). Additionally, the flux 
rope is constructed so that it is not kink unstable. The 
flux rope solution once superimposed onto the coronal 
magnetic fleld leads to an immediate eruption because of 
force imbalance with the ambient magnetic fleld. 

A filament was observed on August 21 (see middle 
panel of Figure 1) and the eruption of the northern sec- 
tion of the fi lament was the cause of the studied eruption 
(|Asai et al.l [2009 ). We place the fiux rope at a position 
and with an orientation in agreement with the observa- 
tion of the filament (see Figure 1). The fiux rope (and 
the line current inside) makes an angle of —15° with the 
—z axis (the Sun rotation axis). It results in a right- 



handed fiux rope with an axial southward magnetic field, 
in agreement with the observations of the filament in Ha 
and w ith the r econstruction of the associated ejecta at 
1 AU ( Asai et al. 2009). The fiux rope is making an an- 
gle of about 20° with the polarity inversion line of AR 
10798 with its negative (resp. positive) footprint in the 
main part of the negative (resp. positive) polarity spot 
of the AR. The direction of the axial field of the fiux 
rope is chosen in agreement with the observat ions of a 
southward (sinistral) filament (Asai et al]l2009l ) which is 
also consistent with the prefilament structure (see Fig- 
ure 4 from iAsai et al.l l2009). The orientation of the axis 
of the fiux rope is consistent with the observations of 
a structure oriented from the northwest to the south- 
east. The exact orientation is chosen so that the positive 
(resp. negative) footprint of the fiux rope are in a region 
of positive (resp. negative) polarity, while the axis of the 
fiux rope remains close to the polarity inversion line. It 
should be noted that in our MHD model, because we use 
a synoptic magnetogram, the positive polarity spot of 
the active region (on the west side of the negative spot) 
is not yet fully developed while the positive polarity on 
the east side is stronger that it is on August 21-22. The 
fiux rope is chosen as right-handed to agree with the 
overall magnetic field in the overlying arcades once the 
direction of the axial field (sinistral) of the fiux rope is 
chosen. This is also consistent with what was reported in 
Asai et al. (2009). Finally, the amount of twist is deter- 
mined by the prescribed toroidal current. The current is 
set by a trial-and-error procedure to match the coronal 
speed of the CME as observed by LASCO (1200 km s-^). 
Because of the way our fiux rope is created (with a dom- 
inant toroidal current), the chosen value of the current 
results in a fiux rope with a much larger amount of twist 
as compa red to what is deriv ed from observations (see 
review by iMackav et alll2010f ). 

Previous time- dependent simulations with the LC 
m odel were performed with the CME initiation model 
of lRoussev et al.l (|2007l ). Here, we use, for the first time, 
a fiux rope model in the LC model because the pres- 
ence of a filament strongly suggests that a fiux rope was 
present prior to the eruption. This CME model has been 
used before to study the evolution of CME in the coron a 
(jLugaz et al.l l2007l : ICohen et al]l2010l : lEvans et al.ll2011f ) . 
A view of the fiux rope at time t = as it is superposed 
onto the steady-state corona is shown with blue field lines 
in the bottom panels of Figure 2. 

2.3. Pre-event Topology 

The magnetic topology before the eruption is very im- 
portant because it determines, to a great extent, the lo- 
cation of the reconnection of the fiux rope with adja- 
cent magnetic field structures during the eruption. The 
steady-state corona around AR 10798 is relatively sim- 
ple with: 1) a large-scale streamer connecting the active 
region with a region of the quiet Sun about 40° east of 
the active region, 2) open field lines with positive po- 
larity, and, 3) active region closed field lines. We found 
no magnetic connectivity between AR 10798 and adja- 
cent active regions, even though AR 10797 is only about 
20-25° away. It is not surprising since AR 10798 is a 
reverse polarity active region (n egative- p ositive) while 
AR 10797 follows the iHale et al.i (jl9T9t ) polarity law 
(positive- negative). Therefore the positive parts of ARs 
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Fig. 1. — Initial configuration of the filament. Left: SoHO/MDI magnetogram. Middle: Ho; image of the filament taken at the Observatoire 
de Paris-Meudon. Right: MHD simulation at time t = with the background magnetic field on the solar surface and the flux rope shown 
with blue fleld lines. In all three views, the images show approximatively 200 arcsec x 200 arcsec. 



10798 and 10797 are next to each other and, moreover, 
they are embedded into positive polarity open magnetic 
fields from the coronal hole. This type of magnetic field 
configuration make it very unlikely that there was a di- 
rect magnetic connectivity between these two active re- 
gions. The topology is much simpler than other cases 
near solar maximum, where Rous sev et al.l ((20071) found 
connectivity between as many as three active regions via 
multiple null points and quasi-separatrix layers. 

There is one null point above the negative polarity 
part of AR 10798 separating the 3 fiux systems as il- 
lustrated in the Figure 2: AR field (pink), streamer field 
(green) and open field (white). Because the negative po- 
larity spot is surrounded by positive polarities, it devel- 
ops into an anemone active regio n (see top right panel 
of Figure 2) a s described before in lAsai et al.l (|2Q08f) and 
iShibata et al.l (|2007|). This null point is originally at a 
height of 0.09 R© above the solar surface. This value is 
probably lower than that on the Sun at the time of the 
eruption because we use a synoptic map of the Sun and 
AR 10798 was not fully developed. The existence of close 
field li nes as part of a unip olar region has been reported 
before (jChertok et al.ll2Q02f ). The steady-state configura- 
tion is an embedded dipole, identical to that of the asym- 
metrical breakout model (Lynch et^ . 2009) and similar 
to that of coronal jet studies (e.g., see lPariat et al.ll2009f ) 
but with a closed outer spine. Here, the inner and outer 
spines are initially closed and there is a simple fan sur- 
face around the negative polarity spot of the active re- 
gion. The inner and outer spines as well as one of the 
lines of the fan surface are shown in red in Figure 2. As 
we superpose the fiux rope onto the active region, the 
null point is pushed by about 4° toward the north-east 
and by 0.04 Rq upward, but the active region retains 
its anemone structure (see bottom panels of Figure 2). 
In addition, a bald patch forms below the fiux rope (see 
bottom panels of Figure 2). Because AR 10798 develops 
and maintains its anemone structure in the simulation as 
was observed on the Sun, we believe our model captures 
the most important features of the solar corona before 
the 2005 August 22 eruptions. 

2.4. Pre- event Comparison of the Simulated Corona 
with EUV Images 

One of the advantages of using the LC model is the 
ability to simulate extreme ultraviolet (EUV) images fol- 
lowing the procedure described in .Downs et al.l ()2010f ). It 



allows us to vahdate our steady-state model by compar- 
ing synthetic images with real ones prior to the event. 
The EUV signal depends on the density and tempera- 
ture of the plasma and it is greatly infiuenced by the 
magnetic structure of the lower corona. Therefore, this 
type of comparisons validates not only the plasma prop- 
erties of our simulated corona but also its detailed mag- 
netic structure. We show a comparison of an 195A and 
17lA images from the EUV Imaging Telescope (EIT, see: 
iDelaboudiniere et al.lll9 95') onboard SOHO 18 hours be- 
fore the eruption with synthetic images for the same fil- 
ters from our simulation in Figure 3. Synthetic and real 
images are plotted using the same scale. The 195A fil- 
ter response function peaks around 1.4 MK, while the 
171A peaks around 1 MK, illustrating different heights in 
the corona. In synthetic and real images, the three active 
regions are clearly visible as regions of enhanced emission 
(hot and dense): the large AR 10800 on the eastern side 
of the Sun near disk center, AR 10797 near the western 
limb and AR 10798 near W50. In addition to the north- 
ern polar coronal hole, there are a number of equatorial 
coronal holes, including two on the south-western side of 
ARs 10798 and 10800. The southern polar coronal hole 
is almost absent in the simulated and real images. Over- 
all, there is good agreement between synthetic and real 
images, which gives us confidence that our model of the 
the corona is a relatively realistic representation of the 
actual corona at the time of the 2005 August 22 erup- 
tion. The most important features for our study that 
the model reproduces are the appearance of AR 10798 
and the presence and aspect of open field regions (dark) 
around it. The main difference between simulated and 
real images is the eastern limb of the Sun where the mod- 
eled emission is too weak as compared to the real one. 
This is relatively unimportant because these regions are 
far from the source region of the CME (more than 120° 
separation) and are not involved in the eruption process. 

3. CME EVOLUTION 

3.1. Early Dynamics of the CME 

As soon as we superpose the fiux rope onto the steady- 
state coronal magnetic field, it erupts due to force im- 
balance. As in other simulations with the same CME 
model, the exact kinematics of the CME early on in the 
corona are not realistic as the CME reaches its maxi- 
mum speed (~ 1500 km s~^) about 1.5 minutes after 
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Fig. 2. — Initial configuration of the corona around AR 10798 from the simulation with the grid structure (top left), the anemone structure 
of the active region before the inclusion of the fiux rope (top right) and after (bottom panels). White, green, pink and blue field lines show 
open, streamer, active region and fiux rope magnetic field, respectively. The solar surface is color-coded with the radial magnetic field 
strength. The bottom right panel shows the zoomed version of the fiux rope as it is added to the solar surface with a view of the bald patch 
(BP) and the null point (NP) as pink isosurfaces. In red, we show the inner and outer spines of the system as well as one of the fieldlines 
of the fan surface. 



the superposition of the flux rope. However, the sim- 
ulated CME kinematics past 3 R© (after 15 minutes) 
are in good agreement with the height-time profile as 
observed in LASC0/C2. Additionally, the CME speed 
after 1 hour is about 1200 km s~^, similar to what is 
measured with LASCO (~ 1250 km s~^). It confirms 
that the total energy of the simulated CME is compara- 
ble to that of the real CME. The model does not intend 
to capture the slow rise phase before the loss of equi- 
librium nor the acceleration phase. This is why when 
comparing synthetic and real images, it is best not to 
use the onset time of the flare as the starting time of the 
numerical simulation but a later time when the CME has 
already significantly accelerated. Previous studies have 
found that the CME acceleration ha ppens during; the X- 
Ray rise phase (Ohvama fc Shibatafi 19981: iFbrb esI 120001: 
iTemmer et all[2008l). We use GOES-12 (jHill et al. 2005) 
and RHESSI CLi n et all l2002[ ) data to investigate the 
flare time in soft X-ray (SXR) and hard X-ray (HXR), 
respectively. For the ejection, the flare onset in SXR was 
00:44UT, the HXR flare started at 01:02UT and the flare 
peaked at 01:22UT. We believe the onset of the HXR flare 
is the best time to use for the start time of our simulation 
since the CME was already observed by LASCO at 4 R© 



10 minutes after the flare peak. 

Figure 4 shows a line-of-sight image of the CME ob- 
served ^^^_JjASC0/C2 and processed with the method 
of Mor gan et al.l ()2006f ) and a synthetic line-of-sight im- 
age from our simulation. The LASCO image is taken at 
01:54 UT, 52 minutes after the onset of the HXR flare, 
while the synthetic image is made 45 minutes after the 
superposition of the flux rope onto the solar surfa ce. The 
synthetic image is processed as explained in Luga z et"al] 
(|2009l) using a synthetic 27-day minimum image created 
from the steady-state simulation. The main difference 
between the synthetic and real images is the latitudinal 
direction of propagation of the CME. In the synthetic 
image, the fastest moving part of the CME is at a PA of 
210° while it is about 220-225° for the real image (this 
region is dimmer than the main part of the CME which 
is closer to PA 270). Because the CME appears as a 
halo, it is not straight-forward to determine its central 
PA but the synthetic image shows a CME propagating 
about 15-20° more towards the south as what appears in 
the real image. The active region was at the Carrington 
latitude South 11 (Sll) when the ejection occurred and 
in our simulation, the CME propagates without a north- 
south deflection, i.e. along a latitude of Sll. As noted in 
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AR 10800 



AR 10798 




AR 10797 



Fig. 3. — Real (left) and synthetic {right) EIT 195A image (top) and 17lA image (bottom) of the corona 18 hours before the start of the 
first 2005 August 22 eruption and from our steady-state simulated corona. On the 171A images, the main active regions (ARs) are pointed 
out in yellow ellipses. 



|Asai et al." (2009) , only the northern part of the filament 
erupted during the flare/ejection studied here, while the 
southern part erupted during the second CME on August 
22. It is therefore likely that the real CME propagated 
more towards the north than what our simulation shows 
since only part of the filament erupted. With our rela- 
tively simple model to start the eruption, it is not yet 
possible to study partial filament eruptions. 

As it rises and expands, the flux rope interacts and 
reconnects with the adjacent magnetic flux systems. Be- 
cause this is a numerical study with finite resolution and 
based on numercial diffusion, the timing of the reconnec- 
tion might not be realistic in our simulation. However, 
we believe that the reconnection process itself reflects 
what occurred on the Sun. It is because: 1) our steady- 
state corona is a good representation of the pre-event 
Sun (see previous section), 2) the CME was due to a fila- 
ment, which was originally at approximatively the same 
position with a similar length, and with the same orien- 
tation as our flux rope, and, 3) the CME kinematics are 
in good agreement with what was observed by LASCO, 
meaning that the total energy in the flux rope is similar 
to that of the CME. In the section below, we analyze 
the interaction of the flux rope with the corona and in 



section 13. 4| we present some observational consequences 
of the reconnection process, which further validate the 
discussed scenario. 

3.2. Reconnection of the Negative Footpoint 

Only a few recent numerical works have focused on 
the interaction of a CME with a realistic magnetic 
structure obta i ned f r om m agneto^ ram measurements 
(iRoussev etla] I20Q7I; iLugaz et al.i 120101 : iCohen et ahl 
2010: ^vans et aLll2011f ). Other simulations have usually 
relied on an ideal representation of the corona with a sin- 
gle dipolar or quadrupolar active refflon in a dipolar Sun 
(as in iManchester 2003; ILynch et al.l 120091 : iJacobs e t al.l 
120091 : IShiota et al. 2010|, for example). Here, the interac- 
tion of the flux rope with the background magnetic field 
is different from both these ideal cases and the realistic 
cases previously studied because of the anemone nature 
of the active region. The topology is different from the 
solar minimum cases studied before with an isolated AR 
not surrounded by equa torial coronal holes (jCohen et al.l 
I2OIOI : lEvans et al.l 120111 ) and also different from the two 
solar maxini u m cas e studies of lRoussev et al.l ()2007l ) and 
"Lugaz et alj ()2010l ) with complex con nectivity between 
multiple ARs. Recently, iTitov et all (|2008. ) studied in 
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Fig. 4. — Line-of-sight image of the CME from LASCO/C2. The left panel is a real image 52 minutes after the start of the HXR flare 
and processed with the normalized radial gradient fllter (NRGF, see: Morgan et al. 2006). The right panel is a simulated image at time 
t = 45 minutes after the superposition of the flux rope onto the solar surface. A simulated 27-day minimum image is subtracted and the 
signal further scaled with distance. 



great detail the magnetic topology of the corona be- 
fore and during the evolution phase leading to the 1997 
May 12 CME. Their study was also for a relatively sim- 
ple magnetic topology without influence from equatorial 
coronal holes. 

The evolution of the CME is altered by the anemone 
nature of the active region, i.e. by the presence of unipo- 
lar magnetic field around it. We find that the positive 
footpoint of the fiux rope does not reconnect away from 
the active region since there is a very limited negative 
polarity magnetic fiux around it. As a consequence, the 
fiux rope remains line-tied at its positive footpoint. In 
contrast, the negative footpoint reconnects quickly and 
extensively with the positive magnetic fiux, partly from 
the positive spot of AR 10798 and partly from the neigh- 
boring coronal hole. The reconnection is a 2-phase pro- 
cess as illustrated in Figure 5. In the first phase (top 
panels), all of the twisted closed field lines of the fiux 
rope (dark blue) reconnect with the streamer magnetic 
field (green). It results in erupting field lines connecting 
the positive footpoint of the fiux rope with the streamer, 
shown in yellow and post-fiare loops in fushia. In the 
second phase (bottom panels), the newly formed (yel- 
low) erupting field lines reconnect with open field lines 
to form (orange) open and twisted erupting field lines. 
Below, we describe the detailed reconnection process and 
the change in magnetic topology during the eruption. 

3.3. Details of the Reconnection Process and Changes 
in the Magnetic Topology 

The reconnection involves two phases: the first one, 
when the negative footpoint of the fiux rope reconnects 
outside the fan surface (top panels of Figure 5) and the 
second one, when some of the erupting field lines open up 
(bottom panels of Figure 5). The detailed reconnection 
process of this first phase goes as follows: first, the (dark 
blue) field lines of the fiux rope reconnect with overlying 
anemone field (pink) to form twisted (light blue) field 
lines connecting the positive foopoint of the fiux rope 
with the negative main spot of the active region (top 
left panel of Figure 5). This reconnection happens at 



the bald patch. In a second step, the newly formed field 
lines (light blue) reconnect with closed (green) field lines 
outside of the fan surface to form the (yellow) field lines 
connecting the positive footpoint of the fiux rope with 
the streamer (top right panel). This reconnection occurs 
at the null point. The two steps of the first phase are 
nearly simultaneous and overlying closed field lines from 
the streamer belt fiux system (green field lines) reconnect 
away. 

During the reconnection process, a null point forms be- 
low the fiux rope and there is a separator linking the bald 
patch and this null point. This separator is similar to the 
bald patch-bald patch separator or the quasi-separatrix 
laye r discussed for exampl e in .Titov fc DemoulinI () 19991 ) 
and lAulanier et al.l (|2QlQf ) and it is shown in metallic 
blue in the middle left panel of Figure 5. A current 
sheet forms in association with the separator and it is 
shown as a purple surface in the bottom left panel of 
Figure 5. The current sheet takes a sigmoidal shape, 
which is very common from fiux em e rgence and shear- 
ing motions (iMaii^ et al.l [2QQ1 iTitov et al.l I2QQ8I : 
lAulanier et al.ll2QlQf ) but it can also be the result of re- 
connection with adjacent fiux systems and due to the 
pre-event topology as is found here. 

During the first phase, which lasts about 5-10 minutes, 
the fiux rope expands to a height of 1 above the solar 
surface. Due to the CME expansion, the null point is 
pushed eastward along the outer spine by about 3° and 
its height increases by about 0.04 R©. At the end of this 
phase, the fiux rope remains composed of closed field 
lines but it appears totally disconnected from its original 
negative footpoint (see top right panel of Figure 6). 

Open field lines (white) from the leading positive spot 
reconnect at the separator below the fiux rope with 
anemone field lines to form new open field lines (shown 
in brown) passing at proximity of the null point. Sep- 
arator reconnection is described in iParnell fc GalsgaardI 
(|2004l ) for example. This is illustrated in the middle left 
panel of Figure 5. A similar type of reconnection fol- 
lowing the eruption of a fiux rope has been previously 
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Fig. 5. — Reconnection process of the flux rope at time 0.5, 1.5, 3.5 (twice), 6 and 45 minutes after the superposition of the flux rope onto 
the solar corona. In aU panels, the sphere is at 1.01 R© and is color-coded with the radial magnetic fleld strength. Dark blue, green, pink 
and white fleld lines show the initial flux rope magnetic fleld, the streamer magnetic fleld, the anemone magnetic fleld and open magnetic 
fleld, respectively. Light blue fleld lines are the result of the flrst step of the flrst phase of the reconnection process. Yellow and orange 
magnetic fleld lines are the result of the 2-phase reconnection process. Yellow fleld lines are closed and connect the positive polarity spot 
of the AR to the streamer belt. Orange fleld lines are open. Fushia fleld lines are post-flare loops. In the middle panels, brown fleld lines 
are newly open. In the bottom right panel, there is an equal number of open and closed twisted fleld lines. See online movies for a different 
view. 



discussed in 'Mac kav fc va n Ballegooiien (2006). Due to 
the reconnection of the anemone flux system and the 
streamer belt and due to the formation of the separa- 
tor, the magnetic topology changes quite drastically from 
the pre-event topology. The null point below the flux 
rope separates the anemone flux system and open field 
lines originating from the leading (white) and trailing 
(brown) positive spots (middle left panel of Figure 5). 
The outer spine of this system is open. Such an opening 
of a similar topology during a time- dependent process 
was shown in lEdmondson et al.l (|2010f ). The eastern null 
point separates the erupting (yellow) magnetic field and 
the closed field originating from the positive trailing spot. 
The outer spine connects the positive trailing spot and 
the positive footpoint of the flux rope through the flux 
rope. The full magnetic topology with all the different 
type of field lines discussed here is shown in the mid- 
dle right panel of Figure 5. The presence of two null 



points separating different but overlapping flux systems, 
as shown here, is only possible in a fully 3-D simulation 
as is the case here. The second phase of reconnection 
happens when the flux rope is not connected anymore to 
the negative spot of AR 10798. 

The second phase of the evolution starts when the neg- 
ative footpoint of the flux rope has fully reconnected out- 
side of the fan surface (light and dark blue field lines no 
longer exist)s. Then, the flux rope has raised, expanded 
and deflected enough so that side reconnection happens 
with open magnetic field lines at the null p oint (similar 
to what was found in lChen fc Shibatall2000l ). The recon- 
nection of closed erupting field fines with open field lines 
from an equatorial coronal hole is a type of interchange 
reconnection. Here, contrary to the case described in 
iCrooker et all (|2002 r). the open field lines are only trans- 
ported over a very small distances as they reconnect with 
closed field lines, because the entire active region is sur- 



SIMULATION OF A CME FROM ANEMONE ACTIVE REGION 10798 



9 




Fig. 6. — Top Left Panel: The pink surface is an isosurface of current density at time t = 5 minutes showing the deformation of the flux 
rope as compared to the initial flux tube shown as a white isosurface of current density. The green fleld hues belong to the streamer belt 
at time t = 0. Top Right Panel: Magnetic fleld lines color-coded with the radial velocity 6 minutes after the superposition of the flux rope. 
This panel illustrates the state of the flux rope after the flrst phase of reconnection. Bottom Left Panel: Initial orientation of the flux 
rope (pink isosurface) with respect to the current sheet (white isosurface). Bottom Right Panel: CME after 1 hour. The pink transparent 
surface is an isosurface of velocity equal to 800 km s"-*^ illustrating the position and orientation of the CME. The white surface is the 
current sheet (same as left panel). It illustrates how the CME is not aligned with the HCS as it leaves the corona. For the left panels, 
the view is approximatively from the ecliptic plane above AR 10798. For the top right panel, the view is from the western limb looking at 
AR 10798. For the bottom right panel, the view is approximatively from Earth. 



rounded by equatorial coronal holes. This special form of 
interchange reconnection is due to the anemone nature 
of the active region and one can expect that it happens 
in most eruptions from an active region inside a coronal 
hole. It has been previously discussed in an observational 
study of an erupt ion from an anen ione active region in 
2007 October 17 (jBaker et al.ll2QQ9[ ). 

At the time the second phase of reconnection starts, 
about six minutes after the CME initiation, the apex of 
the flux rope is at 1 R© above the solar surface and the 
null point is only at a height of 0.2 R©. In this phase 
(bottom panels of Figure 5), the yellow erupting field 
lines connecting the positive footpoint of the fiux rope to 
the streamer reconnect with white open field lines to form 
open, twisted field lines (shown in orange color). During 
this phase, the null point rises and is pushed eastward. 
Its evolution as well as that of the apex of the fiux rope 
is shown in the right panels of Figure 7. Note that at 
all time, the null point is below and to the east side of 
the fiux rope. An online animation shows the first ten 
minutes of the CME evolution with the 2-phase recon- 
nection. 



3.4. Mix of Closed and Open Field and Dimming 
Regions 

The end result of the reconnection described in the pre- 
vious subsection is typically an opening of the fiux rope 
magnetic field but it is a relatively complex process and, 
importantly, some of the magnetic field lines of the fiux 
rope remain closed until the end of our simulation (see 
bottom right panels of Figure 5 and 6). As seen in the 
bottom right panel of Figure 5, at time t = 45 minutes, 
the CME is composed of a mix of closed (yellow) and 
open (orange) field lines. The closed field lines are highly 
twisted and they have, in general, one footpoint (positive 
polarity) near the source region of the CME and the other 
footpoint (negative polarity) in the quiet Sun due to re- 
connection with the streamer belt. The open field lines 
are also highly twisted because they are due to the recon- 
nection of the magnetic field of the fiux rope with that 
of the background originating from the coronal hole just 
south-west of AR 10798 (posi tive polarity). This p icture 
is similar to the sketches of IGosling et al.l ()1995l ) with 
open plasmoid field line embedded inside closed twisted 
field lines. 
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Fig. 7. — Left panel: 5 views of the flux rope at time 1, 2, 5, 10 and 20 minutes after the superposition of the flux rope onto the solar 
surface. Each surface is an isosurface of current density and they are color-coded with the longitude. As the color contour goes from dark 
red to dark blue, the longitude varies from to 10° east of the location of the active region, illustrating the CME deflection. The initial 
longitude of the flux rope in the coordinate system of the Figure is —134.5°; after 20 minutes it is about —143°. This panel also illustrates 
how the CME is deformed but does not rotate. This view is from above the western limb looking down at AR 10798. Right panel: Evolution 
of the apex of the flux rope and that of the pre-existing null point (top: distance; bottom: latitude -red- and longitude -black-). 

gust 24). The fact that dimming from the location of 
the negative footprints of the CME persisted until then 
would corroborate the measurement of a magnetic cloud 
at 1 AU at this time. 



Figure 8 shows a base difference of the 195A image 95 
minutes after the start of the eruption with a pre-event 
image for the simulation (left) and from EIT (middle). 
As highlighted in this Figure, the presence of erupted 
field lines with a footprint at the edge of the streamer 
belt near disk center results in coronal dimmings from 
this region in our synthetic EUV images. Theses coro- 
nal dimmings were clearly visible in EIT images on the 
same region of the solar surface from about 01:26 UT 
on August 22 to the end of August 23. Dimming re- 
gions have been analyzed as the location of the footprint 
of the magnetic cloud (We bb et al.ll2QQQl : iMandrini et al.l 
Hoot) an d/or the region where interchange reconnection 
happens d Attrill et al. 2006, 2007). Here, we find that 
the dimming region corresponds to the footprint of the 
magnetic cloud since erupting closed field lines have their 
negative footprint in this region. We believe that the ap- 
pearance of the dimming region is only possible if the 
first phase of the reconnection process indeed happened 
as described in the previous section. It is hard to under- 
stand otherwise how the filament could reconnect with a 
distant region when reconnection with the adjacent coro- 
nal hole should happen preferentially. 

Our simulation also shows how a fiux rope originat- 
ing from an active region inside a coronal hole can still 
contain closed field lines all the way to the outer corona 
(10 Ro). On 2005 August 24, ACE measured the pas- 
sage of a small magnetic cloud, which ha s been asso- 
ciated with this eruption (jAsai et al.ll2009f ). While the 
smooth rotation only lasted for about 2 hours, it implies 
that twisted field lines associated with one of the two 
eruptions on August 22 made it to 1 AU. Figure 8 also 
includes a base difference image about 7.5 hours after the 
start of the eruption to illustrate the persistence of the 
dimming region. The dimming region was visible until 
the end of August 23 and it is likely that it persisted 
until around 12:00UT on August 24 but the observing 
geometry makes it hard to assess (the dimming region 
had moved close to the western limb of the Sun on Au- 



4. CME DEFLECTION, DEFORMATION AND ROTATION 

4.1. CME Deformation and Rotation 

In the top right panel of Figure 6, we show the fiux rope 
at time t = 5 minutes as a pink surface as well as the ini- 
tial streamer field lines (green) and the initial fiux rope 
position (white surface) corresponding to t = 0. From 
this Figure, it is clear that the fiux rope has been de- 
formed from a straight tube into one with a sigmoid-like 
shape. The reason for the deformation is the reconnec- 
tion with the streamer magnetic field lines (as discussed 
in the previous section). The reconnection happens at 
the negative footpoint of the flux rope (the one on the 
southern part of the sunspot) and it moves the footprints 
of the erupting field westward to the initial position of 
the streamer field lines, as is clear from this Figure. The 
reconnection results in a rotation of the part of the fiux 
rop e which reconnects . It is similar to the rotation found 
in I Cohen et al.l ()2010f ) but, here, we find that it only af- 
fects part of the CME. Also, the net results is not to align 
the CME with the current sheet but to deform it to make 
part of the CME aligned with the streamer belt. Because 
only part of the fiux rope "rotates" , it leads to a deforma- 
tion, which makes the fiux tube looks similar to a sigmoid 
with a marked skewness. This type of skewness has long 
been thought to be associated with shearing motions and 
it is produced by many CME initia tion mec hanisms in- 
volvi ng; shearing: ni otions (iTitov et a l. 2008; Lynch et al.l 
120091 : lAma ri et al"? '2010: Aulanier et al. 2010) and fiux 
emergence (Manchester et al. 2004; Fan 2005). Here, we 
find that a straight flux rope may naturally develop a 
skewness as it interacts with the background magnetic 
field at a height as low as 0.7 R©. This finding is in 
addition to the sigmoidal current sheet developing under 
the fiux rope as discussed in section 3.2. 
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Fig. 8. EUV 195A base difference images of the CME showing the simulated corona at time t = 95 minutes (left) and a EIT image 94 
minutes after the start of the HXR flare (middle) and about 7.5 hours later (right). In all images, the dimming region corresponding to 
the negative footpoint of the closed erupting fleld lines is highlighted by a white ellipse. 



We calculated the orientation of the central axis of the 
flux rope as it propagates from a distance of 1 R© to 
8.5 R0 (flrst hour) and we found no indication of rota- 
tion with an axis making an angle of -15 to -20° with the 
—z axis. The bottom right panel of Figure 6 shows the 
CME after 1 hour as it has reached 8.5 R©. The axis 
of the CME as illustrated by the magnetic fleld lines is 
found to be almost aligned with the initial axis of the 
flux rope. It is also illustrated in the flrst panel of Fig- 
ure 7 with 5 views of the flux rope at different times 
from 1 to 20 minutes (corresponding to distances from 1 
to 3.5 R©). In this Figure, the flux rope is plotted at dif- 
ferent times as isosurfaces of current density. The color 
contours show the CME longitude. The angle between a 
given color and the CME axis gives an indication of the 
CME orientation. The angle is approximatively constant 
over the 20 minutes covered by this Figure, illustrating 
the lack of rotation. T he magnetic fleld m easurements at 
1 AU (see Figure 9 of lAsai et al.l (|2009h ) are consistent 
with a right-handed flux rope with a dominant south- 
ward axial component. It is the same as the orientation 
of the fllament on the solar surface and, therefore, im- 
plies that there was no large-scale rotation during the 
CME propagation. 

Lynch et al^ ([2009|) proposed that the rotation of a 
CME in the corona is related to its skewness. It would 
be due to the force balance required to maintain this 
skewness which is lost as the CME erupts. It generates 
a torque created by the now unbalanced Lorentz force. 
However, here, the skewness is not the result of shear- 
ing motions and it is not associated with force imbalance 
during the CME propagation. Therefore, our conclusion 
that the CME does not rotate is also consistent with this 
scenario. This result provides a caution on the associa- 
tion of skewness with CME rotation as it was found, here, 
that CME skewness can developed during the eruption, 
i n which case n o rotation is expected. 

ICohen et all (|2010) found a strong rotation (- 90°) 
associated with the reconnection of a CME with the 
background magnetic fleld, yielding a CME whose axis 
is aligned with t he heliospheric current sheet (as in 
iShiota et al.ll2010! ). In our case, we found that only part 
of the flux rope reconnected, resulting in a deformation 
instead of a rotation of the flux rope axis. Although re- 
connection continues, it occurs primarily with open mag- 
netic fleld lines, resulting in no change of the CME axis. 



Because there is no observational evidence that the CME 
rotated, it is important to understand why it is the case, 
when most recent simulations have had a tendency to 
show that CME rotation is a very frequent phenomena. 
Here, we flnd that the reconnection of one footprint of 
the CME results in a deformation of the CME but not a 
rotation and that the effect of the magnetic forces is to 
deflect the CME not to rotate it (see next section). 

4.2. CME Deflection 

Figure 7 illustrates the deflection of the flux rope. As 
shown in the right panels, the apex of the flux rope is 
found to deflect by about 10° in 35 minutes as it prop- 
agates to 5 Rq. After this time, the current inside the 
flux rope becomes so small that it is impossible to track 
the flux rope the same way as we did for the flrst 5 R©. 
Instead, we track the position of the velocity disturbance 
associated with the CME. We flnd that the CME further 
deflects by about 1.5° from t = 35 minutes to t = 1 hour, 
at what time the CME has reached 8.3 R©. In the ab- 
sence of a velocity shear in the solar wind which could 
"push" the CME away from the coronal holes, there are 
two possible causes for the deflection: i) the side recon- 
nection with the null point which happens at the east 
side of the CME and ii) the effect of the Lorentz force. 

The flux rope is created by a southward line current 
and the global magnetic fleld above AR 10798 is of unipo- 
lar positive polarity. The Lorentz force acting on the flux 
rope is, therefore, almost purely eastward. The Lorentz 
forc e here is similar to that of the models of lChenI ()l996f ) 
and llsenberg fc ForbesI (|2007[ ) but in these studies, it re- 
sults in acceleration and rotation, respectively due to 
different magnetic topologies. Here, due to the largely 
unipolar polarity, the Lorentz force can deflect the CME. 
Side reconnection resulting in CME deflection can be 
seen in the work of Chen fc Shibata (2000) for exam- 
ple. It is not straight-forward to determine the cause of 
the CME deflection. However, during the flrst phase of 
the reconnection (flrst six minutes), the main part of the 
flux rope appears to be relatively unaffected by the recon- 
nection as its negative footpoint reconnects (see top left 
panel of Figure 7, for example). In addition, the recon- 
nection is a breakout type reconnection (above the flux 
rope) and not a side reconnection. No deflection has yet 
been reported for asymmetric breakout simulations (e.g., 
see .Lvnch et al...2009.) . In addition, the deflection rate of 
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the null point does not correspond to that of the apex 
of the flux rope, as shown in the right panel of Figure 7. 
During the second phase of the reconnection, the side re- 
connection results in an opening of the magnetic field of 
the fiux rope and it is unclear whether it affects or not 
the CME direction. For these reasons, we believe that it 
is likely that the Lorentz force is the main contributor to 
the defiection of the fiux rope. 

The Lorentz force is expected to decrease as the CME 
rises since both the magnetic field strength and the cur- 
rent decrease with height. We find that the deflection 
rate decreases with height and the CME only deflects by 
about 2° in the second 30 minutes of its propagation. It 
is consistent with the Lorentz force being the cause of 
the deflection. The deflection of the actual fllament at 
the Sun could be greater than what is found in the sim- 
ulation because its speed in the low corona, where the 
Lorentz force is stronger, was certainly lower than that 
in our simulation. As noted above, the simulation does 
not reproduce the slow rise of the CME for the 18 min- 
utes from the start of the SXR flare to the start of the 
HXR flare. It might create an additional deflection of 
about 5-10° in the low corona, assuming that the actual 
deflection rate is more or less consistent with what we 
flnd in our simulation. Another reason why our simula- 
tion does not reproduce the observed deflection is also 
due the initiation mechanism. As the force-free fleld of 
the flux rope is superimposed onto the coronal magnetic 
fleld, it results in a magnetic conflguration which is not 
force- free. Therefore, there are additional Lorentz forces 
on all directions, which partially compensate the main 
eastward Lorentz force discussed before. In the real Sun, 
it is likely that the fllament was in a force-free state be- 
fore the eruption and it would only be subject to the 
main eastward Lorentz force, which would yield a larger 
deflection. Deflection of CMEs from the same active re- 
gion during th e following Carrin gton rotation has also 
been reported (jWang et al.ll2006h . 

5. DISCUSSIONS AND CONCLUSION 

We have investigated, using a numerical simulation, 
the flrst of two eruptions on 2005 August 22 from 
anemone active region 10798. We have used a recently 
developed component of the SWMF to include realis- 
tic thermodynamics into the representation of the low 
corona (LC). Using the LC component has enabled us 
to produce EUV images of the solar corona which we 
can compare to FIT observations. In addition, previ- 
ous works have shown that including thermodynamics is 
important to reproduce the plasma and magnetic prop- 
erties of the entire corona. We have investigated the 
pre-event topology of the active region and found that 
it indeed developed into an anemone active region with 
a single null point separating the active region flux sys- 
tem to unipolar positive magnetic fleld. We have found 
that some of these unipolar magnetic fleld lines are open 
while others are part of the streamer belt connected to a 
region near disk center. Pre-event 171 A and 195 A im- 
ages of the corona are in good agreement with simulated 
ones, vahdating our steady-state model. In agreement 
with observations of a Ha fllament, which was well ob- 
served prior to the eruption, we have superposed a highly 
twisted flux rope onto the steady-state solar corona at the 
same location and with the same orientation and chirality 



(right-handed and southward). Due to force imbalance, 
the flux rope immediately erupts and we have focused 
our study on its interaction with the adjacent magnetic 
flux systems and in particular with the open flux from 
the surrounding coronal holes. 

The anemone nature of the source region has two major 
consequences on the evolution of the CME: an eastward 
deflection due to the Lorentz force and a reconnection of 
the negative footpoint of the flux rope, which results in 
the appearance of a long-duration dimming region and a 
mix of open and closed fleld lines in the CME. Since the 
flux rope was southward and right-handed, there was a 
strong southward axial current. The Lorentz force gener- 
ated by this current and the unipolar positive magnetic 
fleld in the vicinity of the active region is east-directed 
and, in our simulation, it deflects the CME by about 
10-15° as the CME propagates to 8 R©. We have not 
address the reason why CMEs from ane mone active re- 
gions tend to be fast CMEs (|Liul l2007i ). It should be 
however noted that, in the simulation, the CME encoun- 
ters closed fleld lines in the low corona, so the magnetic 
topology cannot be fully responsible for the fast speed. 
It is therefore likely that the fast speed of CMEs from 
anemone active regions is at least partially due to the fact 
the CME propagates into a low density , fast solar w ind, 
where the hydro dynamical drag is low ()Cargill|[2004f ). 

The negative footpoint of the flux rope quickly recon- 
nects with the positive fleld part of the streamer belt. 
It leads to a CME with one footpoint (the positive one) 
that is line-tied, and the second one which connects to 
a region of the Sun about 40° away from the active re- 
gion. A dimming region was found to develop in the 
region of the second footpoint, both in simulated EUV 
images and in EIT images and to persist at the Sun for 
as much as a day after the eruption. Additionally, some 
of these closed fleld lines reconnect with open fleld from 
the equatorial coronal hole in a t ype of interc hange re- 
connection previously discussed in lBaker et al.l (|^009). It 
yields a mix of open and closed fleld lines in the erupting 
eject a. While we have conducted our simulation only to 
10 R©, this type of mixed nature of the erupting fleld is 
expected from observations of periods of unidirectional 
strahl inside periods of bi-directional electrons as a mag- 
netic cloud passes oy er a spacecraft (jGosling et al.|[T995l : 
iShodhan et alll2000r ). To the best of our knowledge, this 
is the flrst time such a mixture is obtained from a numer- 
ical simulation of a real CME event. The fact that we 
obtain interchange reconnection in our simulation is not 
surprising, since it has been noted before that anemone 
active regions are the m ost suitable enviro nment for in- 
terchange reconnection (Baker et al.ll2009h . 

We have also found that the CME does not rotate as 
it propagates in the corona (up to a distance of 10 R©). 
This result is found although there is large-scale recon- 
nection, the CME axis is initially not aligned with the 
HCS and the flux rope deforms into a skewed tube at 
a distance of less than 1.7 R©. It shows that the ro- 
tation of CME axis is not a universal process and that 
some CMEs maintain their orientation in spite of large- 
scale reconnection. This result is a pendant to previous 
numerical studies jLvnch et al.l[2QQ9l : iCohen et ani2010l : 
iShiota et al.l l2010f ) which found that rotation happens 
when a CME is initially misaligned with the HCS, or 
when it is skewed. The deflection of the CME, the per- 
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sistence of closed field lines into the inner heliosphere 
and the lack of rotation are also consistent with the ob- 
servation of a short ma gnetic cloud at 1 A U about 2 days 
after the first eruption. lAsai et aTl (j2QQ9f ) associated this 
magnetic cloud with the first of the two eruptions of 2005 
August 22. 

While our simulation captures some of the most impor- 
tant observational features of the CME, there are some 
discrepancies. In our simulation, the CME does appear 
to propagate along the same latitude as that of its source 
region (Sll) while LAS CO images show a CME propa- 
gating close to the ecliptic. One likely explanation is 
that only t he northern part of the filament erupted (as 
reported in lAsai et al.l 120091 ). Our method to obtain an 
eruption is too simple to study partial filament eruption 
or to reproduce the exact dynamics and kinematics of 
the eruption at distances less than 3 R©. Another dif- 
ference between the simulation and the observations is 
the amount of deflection. In our simulation, the CME 
deflects by about 10° eastward whereas a total deflection 
of about 40-50° is expected from observations. A smaller 
deflection might be present in our simulation compared 
to what happened on the real Sun because of different 
kinematics and force balance in the low corona. However, 
we believe that a significant portion of the deflection can- 
not be accounted by magnetic forces, which are dominant 



in the low corona. This finding implies that a significant 
part of the deflection is due to hydrodynamical effect. 
For this CME from an anemone active region, followed 
by a twice faster ejection 18 hours later, these effects may 
take two forms: deflection by the fast wind from the coro - 
nal hole as discussed, for exam ple, in I Wang et al.l ()2004l ) 
and iGopalswamv et aP (|2009f ) and also deflection by the 
second, faster CME pushing the eruption from behind, 
as discussed in IXiong et all (j2009f ). In this case, the ad- 
ditional deflection can be expected to happen more grad- 
ually in the heliosphere. Future simulations will investi- 
gate the heliospheric evolution of CMEs from anemone 
active regions and overtaken by faster eruption. 
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